Mon. Not. R. Astron. Soc. 000, 000-000 (0000) 



Printed 18 November 2011 



(MN style file v2.2) 



O 

> 
O 



00 



Amplitudes of low frequency modes in rotating B type 
stars 



Umin Lee^* 

^Astronomical Institute, Tohoku University, Sendai, Miyagi 980-8578, Japan 



Typeset 18 November 2011; Received / Accepted 



ABSTRACT 

Using weakly non-linear theory of oscillation, we estimate the amplitudes of low fre- 
quency modes in a slowly pulsating B (SPB) star, taking account of the effects of 
rotation on the modes. Applying the formulation by Schenk et al (2002), we compute 
non-linear coupling coefficient between the low frequency modes and estimate the 
i-^ ■ equilibrium amplitudes of the modes excited in the star, assuming the amplitudes of 

the unstable modes are saturated as a result of non-linear coupling with stable modes, 



Q ■ that is, as a result of parametric instability expected between one unstable mode and 

two stable modes. We use the traditional approximation to calculate adiabatic and 
^ ■ non-adiabatic oscillations in a rotating star. We find r-modes in a rapidly rotating 

star play a significant role in the amplitude determination through non-linear cou- 
pling. We also find that for low m modes, the fractional amplitudes of the radiative 
. ^ luminosity caused by the low frequency modes are of order 10~* to 10~^ at the surface. 

^ ! 
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1 INTRODUCTION 

Slowly pulsating B (SPB) stars are a pulsator of low frequency modes excited by the k mechanism associated with the iron 
opacity bump at T ~ 2 x lO^K in the interior (e.g., Dziembowski, Moskalik, Pamyatnykh 1993; Gautschy & Saio 1993). The 
observed amplitudes of the oscillations range from ~ 0.1 to ~ 1 mmag (e.g., Huat et al 2009; Diago et al 2009; Neiner et al 
2009; Cameron et al 2008; Balona et al 2011). Although SPB stars are not necessarily a rapid rotator, the effects of rotation 
on the low frequency modes can be significant, particularly for |2f2/a;| >1, where denotes the rotation frequency and lj is 
d ' the oscillation frequency observed in the co-rotating frame of the stars. 

Here, we are interested in theoretically determining the amplitudes of low frequency modes in SPB stars, taking account 
of the effects of rotation on the modes. Fully hydrodynamical calculation of radial pulsation has a long history, but that 
of non-radial pulsation is not always feasible to carry out, since it usually requires a huge amount of numerical resources, 
particularly for low frequency modes in a rotating star. Instead of hydrodynamical calculation, we may apply weakly non-linear 
theory of oscillation to non-radial pulsation, expecting the amplitudes of oscillation modes are limited by weak non-linear 
coupling between them. As the simplest case, we may consider non-linear coupling between three modes, whose amplitudes 
are expected to reach an equilibrium state as a result of parametric instability between one unstable mode and two stable 
modes (e.g., Dziembowski 1982; Kumar & Goldreich 1989; Wu & Goldreich 2001; see also Craik 1985). 

For the weakly non-linear theory of oscillation, it is essential to calculate the coupling coefficient between oscillation 
modes as well as their excitation and damping rates. To calculate the non-linear coupling coefficient between three oscillation 
modes, we employ the formulation by Schenk et al (2002), who extended the theory to the case of rotating stars. To compute 
oscillation modes in a rotating star, we use the method of calculation by Townsend (2005), who applied the traditional 
approximation to calculate both adiabatic and non-adiabatic oscillations of the star. The traditional approximation is quite 
helpful to largely reduce the computing time spent for calculating a large number of coupling coefficients between various 
non-radial modes in a rotating star. We briefly describe the method of solution in §2, and §3 is for numerical results, and §4 
is for conclusions. 
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2 METHOD OF SOLUTION 
2.1 Linear Oscillation Equation 

Adiabatic oscillation in a uniformly rotating star may be governed by the linear differential equation: 

- + ia;B (0 + C (0 = 0, (1) 

where w is the oscillation frequency in the co-rotating frame, ^ is the displacement vector, B (4) = 217 x ^ with being the 
angular velocity vector of rotation, and C is the differential operation on ^ and its expression as well as its derivation may be 
found, for example, in Schcnk ct al (2002). Assuming the time dependence of the perturbations is given by the factor e"^*, we 
may write the displacement vector as 

i{x, t) = e'^'i(x) = e'"* [C^r + i'ee + ^e^) , (2) 

where Cr, eg, and are the orthonormal base vectors in spherical polar coordinates (r, 9, 4>). We assume that the equilibrium 
of the rotating star is axisymmetric about the rotation axis, and that the dependence of the oscillation modes on the azimuthal 
angle (f> is given by the factor e""'^ with m being an integer representing the azimuthal wave number. The oscillation frequency 
u) may be given as w = a + mfl with a being the oscillation frequency observed in an inertial frame. Because of the dependence 
given by e'(""*+'^*', the oscillation mode with mw < (mw > 0) is a prograde (retrograde) mode. 

Although a separation of variables is in general impossible for oscillations in a rotating star, it becomes possible under 
the traditional approximation, in which the term — sin 9ilee in fl = cos Oiler — sin 6Q,eg is ignored (e.g. , Lee & Saio 1997) . In 
the traditional approximation, the components of ^(a;) are given by 

^=^Wefc^(M;;.)e""^ (3) 

e« = ^^eL(/x;^)e'™^ (4) 
rw^ p[r) 



?* = Z^9lTieL(/";'^)e'"'^ (5) 



P'jr) . 
rw^ p(r) 

and p' and p' , which respectively stand for the Eulerian perturbation of the pressure and the density, are given by 
p' =p'ir)Qkm{t^;u)e'"''^e"^\ (6) 
p' = p'(r-)efc„(yu;i/)e""'^e"^*, (7) 
where p = cos 9, v = 2Q/uj, k is an integer used as a modal index, and 



{l-u^p^)^/l-p^ 



(l-i/V)^!^ 



-vpil - M^)-r- + ^ 
dp 



Qkmip-,!^), (8) 



Qkm{p;y)- (9) 



The function Qkm{p',v), called the Hough function (e.g., Lindzen & Holton 1968), is the eigenfunction, associated with the 

eigenvalue \km, of Laplace tidal equation given by 

Cu [Qkm{p\ I')] = —\kmQkm{p; v) , (10) 

where the definition of the differential operator may be found, e.g., in Lee & Saio (1997). Note that the oscillation modes 
in a rotating star are separated into even modes and odd modes, depending on symmetry of the eigenfunctions about the 
equator of the star. For example, the angular dependence of p' (r,9,4>,t) is symmetric (antisymmetric) about the equator for 
even (odd) modes. In this paper, we normalize the function Q>km as 

/■TT /■27r /•! 

I d9 I d<i> sm.9 \Qkm\^ = 2t^ I dp\@km\' = l, (11) 
Jo Jo J -1 

where 

Bkm = Okme'^'t (12) 

For a given azimuthal wavenumber m, \km depends on the parameter f and tends to lk{lk + 1) with Ik = \m\ + fc as 
— >■ for fc ^ 0, which corresponds to Qkm as — > 0. The quantity V-^fcm represents a kind of surface wave number. 

Except for the prograde sectoral modes (fc = 0), \km increases as v increases. The prograde sectoral modes (associated with 
Aom for modes with moj < 0) are special modes in rapidly rotating stars whose surface wavenumber is lower than the value at 
f2 = 0, and hardly changes with fl. Note that g-modes belong to \km with positive k. On the other hand, r-modes, which are 
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a retrograde mode, belong to X^m. with negative k (Lee & Saio 1997). In the limit of — >■ 0, we have w — >■ 2mfi/i^(i^ + 1), 
where l'). = \m\ + |A; + 1| for negative integer k. Note that Xkm — as w — 2mfl/l'k{l'k + 1)- 

In this paper, wc employ the labeling {Ik, m) with 1^ ~ \m\+k for p-modes and the labeling {1'^, m) with l'), — |m| + |fc+l| for 
r-modes where k is non-negative integer for the former and negative integer for the latter, extending the familiar notation for 
non-radial pulsations of a non-rotating star. Note that in our convention we have mod(/A; — |m|, 2) = and mod(Z^ — |m|, 2) = 1 
for even modes, while mod(Zfc — l^nl, 2) = 1 and mod(Zj^ — \m\, 2) = for odd modes. 



2.2 Weakly Nonlinear Oscillation Equation and Parametric Instability 

Nonlinear evolution of small amplitude oscillation modes in a uniformly rotating star is governed by the oscillation equation 
with nonlinear terms: 

| + B(i)+C(0 = a(^)(tO, (13) 

where ^ = d^/dt and ^ = d^^/dt^, and a^^^ (^)^) represents a collection of nonlinear terms of second order in ^, and the ith 
component of a'^^ is given by (Schenk et al 2002) 

af ^ (1, = -p- V,- {p [(Fi - 1) + Si + mi] } - {l/2)eeV kViVi^, (14) 
where V j denotes the covariant derivative with respect to the coordinate , 

n^ = (v.c^)v.t (15) 

si = (V.$'=)(V,.e^), (16) 

* = (i/2)n[(ri-i)"+ari/ainp] + (1/2) (Fi - i)h, (17) 

n = S'jlli = (V • , H = Sj^i = (Vj^'°) (Vfc^-') , 5j is the Kronecker delta, $ is the gravitational potential, and the repeated 
indices imply the summation over the indices from 1 to 3, and Fi = [dlnp/dln p)^^. Note that we have applied the Cowling 
approximation, neglecting the Eulerian perturbation of the gravitational potential. 

Following Schenk et al (2002), we use eigenvalues uj and eigenfunctions ^ of the linear oscillation equation (1) to expand 
the displacement vector ^{x.t) and its time derivative ^{x,t) in the nonlinear equation (13): 



CA(a3) 

iwA^ytCa:) 



(18) 



for which 

^(CA -iwACA)^A(a;) =0, (19) 

A 

where the subscript A stands for a collection of numbers such as harmonic degree azimuthal order m, and radial order n 
used to identify a linear mode. Note that if a mode with with (ttia, wa) satisfies the linear oscillation equation the mode with 

(— mA, —(jJa) also satisfies the same equation, and both modes arc included in the expansion given above. Substituting the 
expansion (18) into the governing equation (13), and making a scaler product with and integrating over the volume of the 
star, we obtain 

CA{t) - icoACAit) = -i {$A, a'"^ (I, $}) /bA, (20) 
where 

cA{t) = {^A,^Am - im - iB m))) /bA, (21) 

6A = -(lA,iB(lA))+2a;A(lA,^A), (22) 
and for A ^ B we have used a modified type of orthogonality relation given by 

(I A, iB {$b)) - i'^A + wb) {$a, Is) = 0, (23) 
where 



(Ia,Is) = J d''xp{x)eA{x)-$Bix), (24) 

and the asterisk * in the superscript implies the complex conjugation. If we substitute the expansion ^ (x, t) = c% (t)^^ (a;) 
into a'^^ we may obtain 

CA{t) -iwACA{t) = -icJA^ iKABc/eA)c*B{t)ch{t), (25) 

B,C 
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where 

i^ABC = {eA,0.^''HiB,ic)), (26) 

and tA = tOAbA corresponds to the total energy of oscillation in the co-rotating frame of the star (e.g., Lee & Saio 1990). If 
we introduce ca = ca exp(— iwAt), the equation (25) reduces to 

CA = -iwA ^ (nABc/eA) c*B{t)cc{t)e~'^'^\ (27) 

B,C 

where Acu — uja + ^b + i^c- 

If the driving rate of an unstable mode is smaller than the damping rates of stable modes nonlinearly coupled with the 
unstable one, the growth of the unstable mode can be saturated by a transfer of energy to a small number of the damped 
modes. Here, we consider parametric instability between three modes, one unstable mode and two stable modes and we call 
the former the parent mode and the latter the daughter modes, and we expect the amplitude of the parent mode is saturated 
by energy transfer to the daughter modes. In the following, for convenience, we call mode A the parent and modes B and C 
the daughter modes. If we consider nonlinear mode coupling between three modes A, B, and C, we obtain 

CA = -7ACA - il^AVABCCBit)Cc{t)e~'^"* , (28) 

and two similar equations for cb and cc, where we have included the effects of linear dcstabilization (7 < 0) and stabilization 
(7 > 0) of the modes, and we have normalized the eigenfunctions ^g, and such that sa = es = ec = GM^ /R, which 
leads to KABc/eA = kbca/sb = KcAs/ec = r]ABc/'2- 

Parametric instability may occur when the amplitude of the parent mode, |ca|, exceeds the critical amplitude given by 
(e.g., Dziembowski 1982; Arras et al 2003) 



|CA:c|^ = 



\'nABc\'^QBQc 



1 + 



Aw 



7b + 7c 



where Qj = —lOj/'yj. The equilibrium amplitude of the parent mode is then given by 



|CA:e|^ 



(29) 



(30) 



\i1abc\'^QbQc 
and those of the daughter modes are by 

\cB:e\^ ^\cA:ef Qb/Qa, \cc -.e^^ = \cA:e\^ Q C / Q A, (31) 

where A7 = 7a + 7b + 70. Here, we have assumed QbQc > 0, QcQa > 0, and QaQb > 0, which is equivalent to the relation 
given by Qa > 0, Qb > 0, Qc > or by Qa < 0, Qb < 0, Qc < 0, that is, the signs of ujb and uc are the same to each 
other but are different from that of loa, because the parent mode A is assumed unstable (7A < 0) and the daughter modes B 
and C stable (73 > and 7c > 0). Since ujb and ujc have the same sign, to obtain a resonant coupling satisfying Auj ~ 0, 
we have \ujb\ < |wa| and \ujc\ < |wa|. We use the condition A7 > as the criteria for effectively stable equilibrium state of 
three mode coupling (e.g., Wu & Goldreich 2001; Arras et al 2003). 

One of the selection rules giving non-zero coupling coefficient rjABC 7^ is 

niA + TUB + rric = 0, (32) 

and another selection rule may be simply stated that the coupling coefficient ?7abc is non-zero only when the mode triad is 
composed of three even modes or of one even mode and two odd modes (e.g., Schenk et al 2002). For ruA < 0, for example, 
we have two cases because of the selection rule (32), that is, both ms and rnc are positive or one of thb and mc is negative 
so that niBmc < 0. In the former case, if the parent mode is a prograde (retrograde) mode, the two daughter modes are 
prograde (retrograde) modes. In the latter case, however, if the parent mode is a prograde mode having uja > 0, the daughter 
mode with m < is a retrograde mode since ujb < and cue < 0. On the other hand, if the parent mode is a retrograde mode 
having uja < 0, the daughter mode with m < is a prograde mode since ujb > and uJc > 0. 



3 NUMERICAL RESULTS 

As a background model for mode calculation, we use a 4Mq main sequence model computed by a standard stellar evolution 
code, where we have used OPAL opacity (Iglesias & Rogers 1996). The physical parameters of the model are log Tea = 4.142, 
log(L/L0) = 2.470, R/Rq = 2.980, and Xc = 0.4602, and the initial abundance is given by X = 0.7 and Z = 0.02, where T^s, 
L, R, and Xc respectively denote the effective temperature, the surface luminosity, the radius of the model, and the hydrogen 
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mass fraction at the center. Since we use the main sequence model slightly evolved from ZAMS, the model have a thin fi- 
gradient zone above the convective core where /it denotes the mean molecular weight. We employ the method of calculation 
given by Townsend (2005) to calculate in the traditional approximation both adiabatic and non-adiabatic modes of a uniformly 
rotating star, where no effects of rotational deformation are considered. We also employ the Cowling approximation, which is 
good enough for high radial order g'-modes of low degree I. For this model, as well as pulsationally stable low frequency modes, 
we obtain many unstable p-modes and r-modcs, which arc excited by the k mechanism associated with the iron opacity bump 
located at T ~ 2 x lO'^K. It is our main concern here how non-linear three mode coupling determines the amplitudes of the 
low frequency modes in a rotating B-type star. We use the eigenfrequencies and eigenfunctions of adiabatic modes to compute 
the nonlinear coupling coefficient r/ABC- The excitation and damping rates 7 are given by the imaginary part of the complex 
eigenfrequency, uji = Im(a;), which is obtained by non-adiabatic mode calculation. 

In this paper, to prepare a set of daughter modes used to calculate the coupling coefficient rjABC for a given low \m\ parent 
mode, we consider low frequency modes of |m| ranging from \m\ = to 5 and of I in the limited range of \m\ ^ I ^ \m\ + 1 
for |m| and I = 1 and 2 for m = 0, that is, we compute even and odd g-modes with {l,m) = {\m\,m) and (|m| -|- l,m), 
and odd and even r-modes with {l',m) = {\m\,m) and (|m| + l,m) for \m\ = 1 to 5, and even and odd g-modes with 
{I, m) = (2, 0) and (1, 0) for m = 0, in the frequency range 0.05 ^ 2, where ui = ui and M is the mass of the 

star and G is the gravitational constant. Note that the r-modes are in the frequency range of ja;j < 2|m|f2//'(/' -|- 1). For a 
given combination of (mA, mB,mc) for mode triad, the even/odd mode combinations giving non-zero r)ABC are [Ae, Be, Ce), 
{Ae, Bo, Co), {Ao, Be, Co), and {Ao, Bo, Ce), where the subscripts e and o stand for even and odd modes, respectively. 

For a given unstable parent mode, there are numerous combinations of a pair of stable daughter modes, satisfying the 
selection rules for non-zero coupling coefficient tjabc- For each of the combinations we compute the critical amplitude ca.c 
and equilibrium amplitude CA.e for the parent mode using equations (29) and (30), where the excitation and damping rates 
are obtained by non-adiabatic mode calculation. Among the critical amplitudes calculated for various combinations of a pair 
of daughter modes for a given parent mode, we chose the smallest one, considering that the parent mode reaches the smallest 
critical amplitude first to be in an equilibrium state. Since we search for the combination giving the smallest ca-.c from a 
limited set of daughter modes for a parent mode, the critical amplitude ca.c thus determined should be regarded as an upper 
limit for the parent mode. As indicated by equation (29), the critical amplitude \ca:c\ becomes smaller for larger values of 
|??Asc| and y/QsQc and has a dip at Auj — because of the factor 1 + [Auj/{'yB + Jc)]'^ in equation (29). Since normalized 
damping rates 7 = ^GM/R^ of high radial order p-modes of the model are of order 10~^ ~ 10"'*, which are much larger 
than the growth rates ranging from |7| ~ 10~* to ~ 10~® for the g-modes, the dip at Aw = will be very sharp for triads of 
gr-modes having frequencies ^.05. In this paper, we search for the smallest critical amplitude ca-.c for a parent mode among 
mode triads satisfying Aoj ~ 0, that is, nearly in frequency resonance. This procedure is helpful to reduce the number of 
combinations of daughter modes we have to try in order to find the smallest |ca:c|. Because of the assumption Aw ~ 0, the factor 
1-1- [Aw/(7s-|-7c)]'^ in equation (29) takes values between ~ 1 and ~ 10 for most of the triads examined, and \ca:c\ is practically 
dependent on the two quantities |77asc| and ^QbQc, that is, the smallest critical amplitude is likely to take place when either 
|?7asc| or \JQbQg is very large or when both of them are large. Since tiabc, which is calculated by using the eigenfunctions 
of adiabatic modes, is the sum of products of the three eigenfunctions of modes A, B, and C, and since the cigonfunction 
of g-modes has an asymptotic form proportional to cos (J^ krdr) with kr being the wave number in the radial direction 



/(r) being a spatially slowly varying weighting function may be maximized when rig a ~ \ngB — ngc\, which is numerically 
confirmed (see also Wu & Goldreich 2001). The quantity ^/QbQc can be large when a low radial order g-mode or r-mode with 
a very small damping rate is in the mode triad, and in this case we do not necessarily have the property rig a ~ \ngB — Wgc|- 



3.1 Slow Rotation 

In Figure 1, the quantities |r;Asc|, \/QbQc, \cA:e\, and |5I/rad/-f'rad| computed for the mode triad composed of low frequency 
g-modes and corresponding to the smallest \cA:e\ are plotted versus the frequency wa/a/ GM/ of the parent mode for the 



case oiCl = Q,/ ^JGmJr? = 0.01, where 5Z/rad is the Lagrange variation of the radiative luminosity at the surface caused by 

the parent mode, and we have considered only g-modcs for the mode triads, ignoring r-modcs because of the slow rotation. 
Note that we have \cA:e\ ~ |ca:c| for the mode triads plotted in the figure. Since the ratio |2ri/a;| is much smaller than unity 
for unstable g-modes (parent modes) for the slow rotation, the distributions of the points in the figure axe almost symmetric 
between prograde and retrograde modes, although there exists slight deviation from the symmetry. As suggested by the panels 
(a) to (c), \cA:c\ for the parent modes can be small either when |?7asc| or ^/QbQc is very large or when both of them are 
large. For the mode triads in the figure, |77asc| (|cA:e|) tends to increase (decrease) with decreasing |w|, but the dependence 
of |<5Lrad/irad| On u) is not simple because the normalizing amplitudes, defined to make the energy of oscillation equal to 
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Figure 1. Mode triad quantities |r)yiflc'|, VQsQCi |cA;e|i and |(5Lrad/^radl a^r^ plotted versus the frequency u]a_/ \J GM/B? of tlie 
parent mode for a AMq main sequence model for Q./ ^JgmJb? = 0.01, where 5Lrad denotes the Lagrangian variation of the surface 
luminosity I/^ad caused by the parent mode, QbQc = '^B'^c/lBlCy and the red, blue, cyan dots and black open circle stand for the 
parent modes of (/a,™a) = (1, (2, —1), (2, —2), and (3, —2), respectively. Note that for a given parent mode a combination of a pair 
of daughter modes has been chosen so that the critical amplitude |c/i,c| of the parent mode be smallest, and that the thus determined 
combination of daughter modes makes the mode triad associated with the parent mode. Here, the parent modes with positive (negative) 
are prograde (retrograde) modes. 



increase as decreases. The fractional luminosity amplitudes |5-Lrad/irad| at the surface take values ranging from 
~ 10"" to ~ 10~*, and those of the parent Z = |m| = 1 g-modes tend to be smaller than the others. 

As shown by Figure 1, the coupling coefficient \rjABc\ can be as large as ~ 10® for the very low frequency parent modes 
for n = 0.01. In Figure 2, the coupling coefficient rjABcij') — kabc{i') / (see Appendix B) and the eigenfunctions xzi and 
XZ2/ iciu)'^) are plotted versus x = r/R for a mode triad composed of a parent mode of Ia ~ —rriA = 1 and daughter modes of 
Ib ~ —ruB ~ 4 and Ic ~ mc = 5 for Q, — 0.01, where the eigenfunctions are normalized so that the oscillation energy be equal 
to GM^ /R. The panel (a) shows that the mean magnitude of IriABcir)] increases rapidly with increasing r in the ^-gradient 
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Figure 2. Coupling coefficient »7asc{'') = f^ABci^) / and eigenfunctions xzi and xzij {c\G)^') versus x = r/R for a mode triad 
composed of a parent mode of lj\ = — = 1 and daughter modes of Ig = —nig = 4 and Iq = mc = 5 for fj = 0.01, where the 
eigenfunctions zi and Z2 are normalized so that the oscillation energy in the co-rotating frame be equal to GM^ / R with G, M and R 
being the gravitational constant, and the mass and radius of the star, respectively. See the Appendix B for the definition of the quantities 
ci, zi and Z2. Here, Gjj\ = 0.1680, uib = —0.07082, and ujc = —0.09769. In panels (b) and (c), the black, red, and blue lines indicate the 
mode A, B, and C, respectively. 



region above the convective core, and the panel (c) indicates that this rapid increase is caused by the amplitude trapping of the 
eigenfunctions into the /i-gradient zone. Note that rjABcir) stays almost constant outside the /.i-gradient region as shown by 
panel (a). For this mode triad, the modes B and C are very high radial order jr-modes, for which we have UgA ~ \ngB — ngc\- 
The panel (a) suggests that for the mode triad the non-linear coupling between the modes preferentially occurs in the thin 
/x-gradient zone having a high Brunt- Vaisala frequency. Note that for the ZAMS model with no /^-gradient zone, there occurs 
no rapid increase in the mean magnitude of \'riABc(T)\ immediately above the convective core. 



3.2 Rapid Rotation 

For weakly non-linear coupling of oscillations in a rapidly rotating star, r-modes may come into play as an important member 
in mode triads for parametric instability. For a rotation rate Q. = 0.2, for example, the maximum frequency 2\m\fl/l' {I' + 1) 
for even (odd) r-modes in the co-rotating frame is 1/15 (0.2), 1/15 (2/15), 0.06 (0.1), 4/75 (0.08), 1/21 (1/15) for \m\ = 1 
to |m| = 5, and the corresponding inertia! frame frequency \a\ is 2/15 (0), 1/3 (4/15), 0.54 (0.5), 56/75 (0.72), and 20/21 
(14/15), respectively. In Figure 3, the mode triad quantities |77asc|, VQbQc, \cA:e\, and |(5Lrad/i'rad| are plotted versus the 
frequency uja of the low \m\ parent modes for fl = 0.2. The distribution of the points of the parent modes of I = \m\ = 1, for 
example, is not symmetric any more between prograde and retrograde modes, because the frequency spectra of low frequency 
(;-modes themselves largely deviate from the symmetry for |2n/cj| ^1, and r-modes appear only as a retrograde mode. Since 
the damping rates 7 of low radial order r-modes can be smaller than lO"^'' and those of low radial order 3- modes are of order 
10~* to 10~^, it is likely that the mode triads giving the smallest critical amplitude \ca:c\ for the parent modes are those 
containing a low radial order r-mode or jr-mode. Figure 4 shows mode triad quantities |r;Asc|, VQbQc, and |5irad/irad| the 
same as those plotted in Figure 3, but here the red, blue, and cyan dots respectively stand for the mode triads containing 
no r-modes, one r-mode, and two r-modes. The low frequency parent modes tend to be coupled to one r-mode or two. We 
find cases where the parent modes are an unstable r-mode coupled with a stable r-mode. Since both |7?asc| and ^/QbQc are 
large, the amplitude |(5Lrad/irad| becomes very small. We also find cases in which the parent modes are a low frequency even 
g-mode coupled with two stable odd r-modes, one low radial order and the other high radial order r-modes. If no r-modes 
are in a mode triad, on the other hand, the amplitude I^Lrad/^radl for the parent mode is in general larger than those of the 
mode triads that contain one r-mode or two. 

In Figure 5, the quantities rjABcir), xz\, and xz2/{c\ijj^) are plotted versus x — r/R for a mode triad composed of a 
parent mode of Ia ~ —rriA ~ 2 and daughter modes of /s — 1 = ms — 1 and Ic — 1 = rnc = 1 for H = 0.2, where the 
eigenfunctions are normalized so that the oscillation energy be equal to GM^ /R. Here, the daughter mode C is an r-mode, 
and the eigenfunction zi of the r-mode has an amplitude much smaller than those of the (j-modes A and B, although the 
amplitude Z2/{c\u]^) of the r-mode is comparable to the g-modes. The figure suggests that the terms containing Z2/{c.\uj^) 
and/or d[z2 / {ciCiP)]/ dx make dominating contributions to tjabc (see Appendix B). We also note that the amplitudes of the 
modes in the triad are not strongly trapped in the /i-gradient zone. Although the coefficient riABcix) is spatially oscillatory 
in the region 0.1 ^a; ^0.3, it takes almost a constant value in x ^.5. 
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Figure 3. Same as Figure 1 but for the case of r2 = 0.2. Here, we have included r-modes as a possible member in the mode triads to 
calculate the coupling coefficient »7asc- 



3.3 In an Inertial Frame 

It may be useful to plot the fractional amplitude of both the parent and daughter modes as a function of the oscillation 
frequency observed in an inertial frame, where the daughter modes are regarded as being non-linearly excited by the parametric 
instability. Figure 6 shows |5I/rad/^rad| versus the inertial frame oscillation frequency for the case of H = 0.01 (panel a) 
and n = 0.2 (panel b), where a — cj — mCl, and the filled (open) symbols stand for the parent (daughter) modes. Here, only the 
daughter modes having Z ^ 3 are plotted in the figure. If one of the daughter modes in a mode triad has a damping rate much 
smaller that the excitation rate of the parent mode such that \Qb\ 3> \Qa\ or \Qc\ 3> \Qa\ (see equation (31)), the amplitude 
|5^rad/irad| of the daughter mode can be comparable to or even larger than that of the parent mode. In fact, although the 
upper limit of |5Lrad/irad| for the parent modes is of order ~ 10"**, there are many daughter modes whose amplitude is as 
large as |5I/rad/^radl >10~'^. Since \ujb\ + |iic| <l'^.4|, the daughter modes are likely to be in the low frequency domain of 
\a\ for the case of f2 = 0.01 since the term mQ, in a — ui — mfl is small compared to u and hence tr ~ a). For the case of 
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Figure 4. Mode triad quantities |r;ABcl i VQbQc^ a^nd I "S^rad /^rad I a-re plotted versus the frequency uja/ ^JgM/B? of tfie parent mode 
for a 4M0 main sequence model for the case of f2 = 0.2, where (SLj-ad denotes the Lagrangian variation of the surface luminosity Lj-ad 
caused by the parent mode, QbQc = '^B'^c/lBlC^ and the red, blue, cyan dots indicate the mode triads containing no r-modes, one 
r-mode, and two r-modes, respectively. Here, the parent modes with positive (negative) are prograde (retrograde) modes. 
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Figure 5. Coupling coefficient riABc(^) = f^ABc(^)/^A and eigenfunctions xzi and XZ2/ (ciuP) versus x = r/R for a mode triad 
composed of a parent mode of 1^ = —ttla = 2 and daughter modes ol Ig — 1 = mg = 1 and Iq — 1 = mc = 1 for f2 = 0.2, where the 
eigenfunctions 21 and 22 are normalized so that the oscillation energy in the co-rotating frame be equal to GM^ /R with G, M and 
R being the gravitational constant, and the mass and radius of the star. See the Appendix B for the definition of the quantities ci, z\ 
and 22. Here, u)a = —0.5332, u)b = 0.3329, and Coc = 0.1999, and the mode C is an r-mode. In panels (b) and (c), the black, red, and 
blue lines indicate the mode A, B, and C, respectively. 



n = 0.2, however, the term mCl can be comparable to Oj and even the daughter modes can have the inertial frame oscillation 
frequencies comparable to those of the parent modes. The range of the fractional amplitude |(5Lrad/irad| for Q, = 0.2 is much 
wider than that for Q, — 0.01, although the upper limits are almost the same. It may be interesting to note that for the 
case of n = 0.2, low radial order odd I' = m = 1 r-modes, which are linearly stable but non-linearly excited, have very low 
oscillation frequency o" ~ in the inertial frame, although the amplitudes are not very high. As indicated by the existence of 
vertical sequences of open symbols (daughter modes) in the panels (a) and (b), there arise some cases in which a stable low 
radial order r-mode or g-mode, which has a very small damping rate 7, is shared by several parent modes, indicating that one 
stable daughter mode has different amplitudes depending on the mode triads it belongs to. This may suggest that the set of 
daughter modes we use for the computation of rjABC is not large enough, or that the three mode non-linear coupling theory is 
too simplified to be applied to the case where modes having an extremely small damping rate exist in dense frequency spectra 
of oscillation modes. 



4 CONCLUSIONS 

Using the weakly non-linear theory of oscillation, we have estimated the amplitudes of low m y-modes and r-modes in a slightly 
evolved AMq main sequence model, assuming the mode amplitudes are limited by parametric instability between one unstable 
mode and two stable modes. Here, the unstable low frequency modes are assumed destabilized by the K-mechanism associated 
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Figure 6. Fractional amplitude |<5irad/^radl of the surface luminosity versus the oscillation frequency \a\/ ^yGM/I{^ observed in an 
inertia] frame for f2 = 0.01 in panel (a) and for S7 = 0.2 in panel (b), where the filled (open) symbols stand for the parent (daughter) 
modes, and the legend of the symbols for both panels are given in the lower right corner in panel (a), and only the modes of 1 ^ i 3 
are plotted. For the parent modes, the combination {I, |m|) is given by (1, 1), (2, 1), (2, 2), and (3, 2), while for the daughter modes the 
combination (Z, |m|) is given by (1,0), (1,1), (2,0), (2,1), (2, 2), (3,2), and (3,3). 



with the iron opacity bump, and we have taken account of the effects of rotation on low frequency modes in the traditional 
approximation. For a given unstable mode (parent mode A), we compute three mode non-linear coupling coefficient rjABC 
for various combinations of two stable modes (daughter modes B and C), and we choose, among the numerous combinations, 
the one that gives the smallest critical amplitude. The traditional approximation is employed in order to reduce the amount 
of computing time necessary to find the optimal combination of daughter modes. It is important to note that since we can 
use only a limited set of daughter modes the critical amplitude thus determined for a parent mode should be regarded as 
an upper limit. The critical amplitude essentially depends on |?7asc| and \/QbQc if we assume resonant mode coupling 
satisfying Au ~ 0, and the smallest critical amplitude may take place when either |r;Asc| or ^/QeQc is very large or when 
both of them are large. If the damping rate of a parametrically excited daughter mode in a mode triad is less than the growth 
rate of the parent mode, the equilibrium amplitude of the daughter mode can be larger than the parent mode. It is therefore 
likely that parametrically destabilized daughter modes like low radial order jr-modes and r-modes are among the periodicities 
observed in a rapidly rotating B star. The fractional amplitudes 1 5Lrad /-f'rad | of the parent and daughter modes can be of 
order ~ 10"* to ~ 10"^ for the main sequence model, the magnitudes of which may be consistent with those observed in B 
type variable stars (e.g., Huat et al 2009; Diago et al 2009; Neiner et al 2009; Cameron et al 2008; Balona et al 2011). We 
also find that the amplitudes j^-Lrad/iradj of the parent mode tend to be large for high I values, although the visibility of the 
modes decreases with increasing I. 

We find that r-modes significantly affects the amplitude determination of low frequency modes in a rapidly rotating star. 
Since low radial order r-modes of the model have damping rates -y as small as or even smaller than ~ 10"^", the mode triads 
giving the smallest critical amplitude for the low frequency parent modes are likely to have a low radial order r-mode as 
a member. When a mode triad has a low radial order r-mode, the equilibrium amplitude of the parent mode tends to be 
smaller than those for the mode triads without r-modes. We find some cases in which a low radial order r-mode is shared 
by several parent modes. We think the degeneracy of the daughter mode in mode triads is a problem in the weak non-linear 
coupling theory we use, since the equilibrium amplitude of the daughter mode depends on the parent modes it is coupled 
to. This degeneracy might be removed if we use an enlarged set of daughter modes to determine the optimal combination, 
or this degeneracy may suggest that the three mode coupling theory we use is too simplified to be applied to the problem 
we are considering, that is, we have to consider higher order non-linear mode coupling to lift the degeneracy. In spite of the 
problem, the weakly non-linear theory could be useful when we try to compare theoretical mode calculations to observations. 
Generally, the number of observationally detected low frequency modes for SPB stars is much smaller than that of theoretically 
calculated linearly unstable modes (e.g., Walker et al 2005, Saio et al 2007), and we may use the weakly non-linear theory to 
decide which linearly unstable modes can have amplitudes large enough to be detected observationally. We may suggest that, 
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if a low radial order I' = m = 1 r-mode, which may be linearly stable, is parametrically excited, the r-mode can produce 
very long period variations observed in an inertial frame. We may attribute very slow pulsations detected in SPBe stars (e.g., 
Walker ct al 2005, Saio et al 2007) to low radial order I' = m = 1 r-modes, although the amplitudes would not be very high 
as indicated by Figure 6. 

We have carried out an additional calculation to obtain the optimal critical amplitudes \ca:c\ for the parent modes by 

extending the set of daughter modes from /max ~ \m\ + 1 = 6 to /max = 8 for H = 0.01, and wc obtained almost the same result 
for their amplitudes |(5Lrad/irad|. However, it is extremely time consuming to carry out similar calculations for /max much 
larger than /max ~ 10, and from the numerical results we currently have it would be fair to say we are not able to correctly 
specify what the most likely degrees of the daughter modes are for the parent modes. To construct a set of daughter modes 
extended for a very large /max, asymptotic methods would be useful to represent the eigenfunctions and eigenfrequencies 
for high / and high radial order g-modes and to calculate the coupling coefficient tjabc (e.g., Dziembowski 1982), and we 
may use an asymptotic treatment by Lee & Saio (1989) for low frequency modes in uniformly rotating stars. Extending our 
weakly non-linear analysis to large values of /, we can also consider non-linear couplings not only between a parent mode and 
many pairs of daughter modes having similar frequencies but also between a daughter mode and granddaughter modes with 
frequencies still lower than that of the daughter mode (e.g., Kumar & Goodman 1996). The problems of these highly multiple 
mode couplings could be important for the amplitude determination for both the parent modes and daughter modes, and we 
may have to include these mechanisms in our analysis to obtain definite answers for the oscillation amplitudes in the stars. 

We have assumed uniform rotation in the present analysis. It is, however, quite likely that differential rotation is a rule in 
reality in a rotating star, and that even a weak differential rotation would affect the frequency spectrum of the low frequency 
modes. We need to understand the property of low frequency modes in a differentially rotating star as well as how a differential 
rotation law is established in a star. The problem of differential rotation in a rotating star is quite difficult to find answer and 
is beyond the scope of this paper. Since we used the traditional approximation to compute low frequency modes in a rotating 
star, we could not correctly take account of the effects of linear coupling between the modes associated with different \kmS- As 
discussed by Aprilia ct al (2011), the linear coupling between low frequency modes tends to preferentially stabilize retrograde 
g-modes for rapidly rotating B stars, particularly for those having lower effective temperatures. It is therefore desirable to 
use the expansion method (e.g., Lee & Saio 1987; Lee & Baraffe 1995) to compute low frequency modes in a rotating star for 
the weakly non- linear coxipling calculation (as well as for the analysis of the pulsational stability), although the calculation 
using the expansion method would be much more time consuming than that using the traditional approximation to find the 
optimal combination of daughter modes for a given parent mode. In a B type main sequence star, a jtt-gradient zone with a 
high Brunt Vaisala frequency forms above the convective core as it evolves from the ZAMS, and the /i-gradient zone has the 
effect of enhancing the coupling coefficient |»7abc| compared to the case of the ZAMS model with no /i-gradient zone. This 
enhancement of |?7asc| would affect the equilibrium amplitudes \cA:e\ of the low frequency modes. It is therefore important 
to examine the effects of stellar evolution on the quantities |r/ABc| and |cA:e| and hence on the amplitude determination of 
the low frequency modes in SPB stars in the weakly non-linear coupling theory. 



APPENDIX A: OSCILLATION EQUATION IN ROTATING STARS IN THE TRADITIONAL 
APPROXIMATION 

The oscillation equations for uniformly rotating stars in the traditional approximation may be given by (e.g., Lee & Saio 

1990) 

= (ciu'' + rA)zi + il-U-rA)z2, (A2) 
dr ^ ' 

where 

., = m, = (A3) 

r pgr 

and 

^^_dlnp^ dlnM. ^ ^ cilnp _ ^ dlnp 
dlnr' dlnr ' dlnr Fidlnr' 

ci = {r/Rf'/{Mr/M), Mr = J^iTrr^pdr, g = GMr/r^, Ti = (dlnp/dlnp)^, and M and R are the mass and radius of 

the star, and G is the gravitational constant. With appropriate boundary conditions imposed at the centre and surface of 
the star, we solve the above set of differential equations as a boundary-eigenvalue problem for the eigenfrequency oj. Since 
^'^ = rziQkm, = ^-?26fem/cia;^, and ^'^ = rz2B'^^/ciui^, the derivatives d^J" /dr, d^^ /dr, and d^^/dr may be calculated by 
making use of equations (Al) and (A2). 
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Using the eigenfunctions ^, the oscillation energy e observed in the corotating frame of the star is given by (Lee & Saio 
1990) 

rR 

e = a)6 = / ^* • ^pr^dr, (A5) 



and it is interesting to note that for positive \kn 



I = w / 
Jo 

is interes 

(a,^ + ^^) fdt,\Qkm\'' = fd„[\Ql4+HJ\)- (A6) 
APPENDIX B: CALCULATION OF COUPLING COEFFICIENT kabc 

If we neglect the boundary terms using the pressure zero surface boundary condition, by use of partial integrations we can 
rewrite the expression for the coupling coefficient kabc as (Schenk et al 2002) 



where 

.(1) - i / ^3 



^ABC 



2 



I^ABC = o / d^XP 



d'xp (Fi - 1) (Has V • + "^bcV ■ ^a + ^caV ■ ^g) , (B2) 
V • ^aV • IsV • ^c, (B3) 



(ri-i)^+ 



d\np 

ABC = \ j d^^P iXABC + Xacb) , (B4) 

,(4) _ _1 / ,fi^.-.c\ci ffc' 



^ABc = -2 y d'^a^pCAfBC&V.V.Vfe*, (B5) 

where Eab = <5]S^' (^a,^s), Xabc = S^Xj (^a,4b,4c)> (Ia,Ib,Ic) = (ViCA)(VfcCs)(ViCc), and the repeated indices 
imply the summation over the indices from 1 to 3. We note that since xabc = Xbca = Xcab and xacb = Xbac = 
XCBA, the coefficient kabc does not depend on the order of the indices A, B, and C. Note that if we consider kasc(^) = 
(dKABC /dr)dr, we have kabc = kabc{R)- 
In spherical polar coordinates (r, 6, 4>), the covariant derivatives of the displacement vectors are 

«r-f-^*. (Be, 

e.-f -f-(«^e-^6»), (BT) 

e^ = ^- sin = ir (z,mQ - ^ sin ee^) , (B8) 
d(p \ cioj^ J 

^■,r = = -^[r—-AQ , (B9) 



r dr r dr \ cilo'^ 

= + ^- = z,Q+^^dQ\ (BIO) 

r ad r ciuj-' 

v<A = -^ cos6'^'^ = i — r^B^, (Bll) 



rsmo or or \ ciu>-' J smt/ 



r sin 6 86 



e,^-L-% + il + i£2£|^« = + ^e^, (bi4) 

•® rsmy ocp r rsm^ ciuj^ ^ 



where we have written ^'j, instead of VjC, for the covariant derivatives, and 
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and 

ei = me*-cos0e^ el = -JI!-e^ + S^q'>. (bi6) 

sin y sin 9 

Note that wo have omitted the subscript km attached to the functions 0km, ©fcmi ^tm' ^^^-j for simplicity. 
The integrands of k may be given by 

e j + 5 {ziziH : eee) + s {j^^iH ■. [ae*G + e^e] (bit) 



2 \ ClLO^ ClUl^ 



sin 6 



Y \or ciLO'' J sm6 J \ ciu)'' \or c\u-^ ) ^ 

V • IaV • ^bV • ic = HaHbHcQaQbQc, (B18) 

+s ( ^„ (#r^) : e faee* - _ s [e'e- - e'e'l ) 

y or \dr c\u^ J \^ smd J ^ ^ j 

+5 ( ) : dQdQ'Q" ^ ^ ^ e e*e' - e-^e-^ 

I c\u^ \or ciuj-' J smd 

KciLj-' cioj-' \or ciw^/ ^ sine/ y 

_.c/ /"^ «2 \ ;,/snne ^eee-^^i (B19) 

+S zizi — — : 9666 - ' 



sm d 

I ciw^ ciw^ ciw^ ' sin^ I \^ ciw^ ciw^ ' sin^ 

+^5 fzi^^ : ee^e^ + eae'se") + \s (zizi^ ■. eeae* + eee^) 

/ d{rzi)A d{rzi)B d{rzi)c , ^ . ^ , ^ \ A r> A 
\ — dr 2(zi)a{zi)b(zi)c I 6.6360 



+ 
and 



r-^^kkfe*... = Is (-1^^ : e [e^e« - e^e^]) |: (^f ) + (.i).(.i)B(^i)c6.eBec||, (B2o) 

where the function H{r) is defined by 

V ■ 4 = i/(r)e(0, ^) = -^{z2- zi) e{9, 0), (B21) 
i 1 

and 

,1 J.2 J.3 12 3\ /•! r2 r3 1 2 3 , rl r2 r'i 1 2 3 . rl J.2 r3 1 2 3 



CY / /•! /•2 J.3 1 2 3\ J.1 r2 r3 1 2 3 , /•! (-2 r'.i 1 2 3 , rl /■2 r3 1 2 3 

S [J J J -P P P ) = JaJbJcPaPbPc + JaJcIbPaPcPb + JbJcJaPbPcPa 



I /■l/"2r3 1 2 3 1 p\ c'2 p'A 1 2 3 i /■l/'2r3 1 2 3 

+JbJaJcPbPaPc + JcJaJbPcPaPb + IcJbJaPcPbPa, 



s {fff : pW) dn = fAflfcZ'^ic + Afcflz'iSB + fifcfAZ^icA 

. fl /2 j3 7I23 , fl f1 j-3 ^123 , i-1 rl c3 7-123 

'tJbjajc-^bac + JcJaJb^cab + JcJbJa^cba-, 
where 



(B22) 



where the functions /-'s depend only on r and the functions p-'s only on S and (j). Note that 

S (/VV^^ : P W) = S (/W^^ : pVp') = S [fff : pW) = 5 (/VV^ : p'pV) = ■ ■■ (B23) 
Integrating S [f^ f^f^ ■ p^p^p^) over a sphere of radius r, we obtain 



(B24) 



Zabc = / PAPBPcd^. (B25) 
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